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Abstract. The non-vanishing of naive T-odd parton distributions function can be explained by the existence of the gauge 
link which emerges from the factorized description of the deep inelastic scattering cross section into perturbatively calculable 
and non-perturbative factors. This path ordered exponential describes initial / final-state interactions of the active parton 
due to soft gluon exchanges with the target remnants. Although these interactions are non-perturbative, studies of final state 
interactions have been approximated by perturbative one-gluon exchange in Abelian models. We include higher-order gluonic 
contributions from the gauge link by applying non-perturbative eikonal methods, incorporating color degrees of freedom in 
a calculation of the Sivers function. In this context we study the effects of color by considering the FSIs with Abelian and 
non-Abelian gluon interactions. We confirm the large Nc QCD scaling behavior of Sivers functions and further uncover the 
deviations for finite A'e. Within this framework of FSIs we perform a quantitative check of approximate relations between 
T-odd TMDs and GPD which goes beyond the discussion of overall signs. 
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Over the past two decades the transverse partonic structure of hadrons has been the subject of a great deal of 
theoretical and experimental study. Central to these investigations are the observations of large transverse single spin 
asymmetries (TSSAs) and azimuthal asymmetries in hadronic reactions- from inclusive hadron production [1, 2, 3] 
to Drell-Yan Scattering [4, 5], and in semi-inclusive deep inelastic lepton-nucleon scattering experiments [6, 7]. Two 
explanations to account for TSSAs in QCD have emerged which are based on the twist-three [8, 9, 10] and twist- 
two [11, 12, 13] approaches. We focus on the twist-two approach in the factorized picture of semi-inclusive deep 
inelastic lepton-hadron scattering (SIDIS) [13, 14] at small transverse momenta of the produced hadron, Pt ^ kj « 
\f^, where is the hard scale. In this kinematic regime the Sivers effect describes a twist-two transverse target 
spin asymmetry through the "naive" time reversal odd (T-odd) structure, A/(ji[:,^7') ~ 5r • x kT)fij{x^k^) [11]. 
kj- is the quark intrinsic transverse momentum and P is the momentum of the target. The Sivers asymmetry has 
been the focus of much theoretical work on QCD factorization theorems. Among the most interesting results is 
that the Sivers function is not universal. It is predicted that there is a relative sign between the Sivers function from 
consideration of the gauge link dependence going from SIDIS to Drell-Yan scattering [15, 16]. The non-universality 
is further reflected in azimuthal asymmetries where one encounters the transverse moments of the quark correlator 
<I> [17]. Here <t>'^(x) — J d^kj kj<P{x,kT), where the non-trivial link dependence remains after integration over kj- 
The correlator decomposes as <I>" (x) = <i>" (x) + Cg 7Z'Pq{x,x) with calculable process-dependent gluonic pole factors 

Cg [18] {e.g.Cc-^"^'^^ = -Cg*^^-' whereas Cg^^"^'^^ = !)■ contains the T-even operator combination, while 4>g 
contains the T-odd operator combination. The latter is precisely the soft limit xi — >^ of a quark-gluon coiTelator 
<I>g('^7'^i)- It was shown in Ref. [17] that the fej-weighted Sivers function can be written in terms of the gluonic pole 
matrix element which we express 

(fc7>W = / |^e«^"^7P,5r|^(-z-n/2)7+[-z-«/2;z-«/2]/'(z-«/2)^(z-«/2)|P,5r), (1) 

where {kj){x) = 2Mei^s\.f^j^\x). The light-like vector n = (n^ = l,n+ = 0,«x = 0) represents a specific direction on 
the light-cone, and [x ; y] denotes a gauge link operator connecting the two locations x and y. The operator / originates 
from the time-reversal behavior of the FSI/ISI implemented by the gauge link operator in (1) and is expressed in terms 
of the gluonic field strength tensor 



2p{z-n/2) = / dy-[z-n/2,y-n]gF+'{y-n)[y-n,z-n/2]. 



(2) 
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FIGURE 1. Left: The matrix element W = {P — k\['x> n;0]q{0)\P) dressed with the FSIs. The FSIs are described by a non- 
perturbative scattering amplitude M that is calculated in a generalized ladder approximation [20]. Right: The quark-quark correlator 
with FSIs. 



We exploit the model dependent relation of transverse distortion in momentum space with transverse distortion in 
impact parameter space [19] by inserting a complete set of momentum states and demand that the operator / is diagonal 
in momentum eigenstates, that is La^Lap • • • (■^p|^'|'^p)^'^H'^p ^ ■^p)- This leads us to model the target remnant of 
the quark-quark correlation function in terms of sum of spectators states {|/lp)} and to define the lensing function 
/(P) = (■^p|/|'^p)5'''^'' (Ap — Ap). In this picture is it clear why spectator models can describe the gluonic pole matrix 
element as the factorization of FSIs convoluted with transverse impact parameter distortion (via impact parameter 
GPDs). While one cannot achieve a general factorization between TMDs and impact parameter GPDs; that is for the 
whole sum, the phenomenological relation may exist for each Fock space state separately. In this approximation 
the Dirac structure is carried by the distortion and the color dependence is determined by the chromodynamic 
lensing function [21]. This is most clearly displayed by transforming Eq. (1) into a mixed coordinate-momentum 
representation where the gluonic pole matrix element is expressed as 

(k-rKx) = I d^br I ^e«^^^"(P+,0r;5r|^(zi)7+[zi;z2]/'fe)?fe)|^+,0r;5r). (3) 

The impact parameter iij- sits in the arguments of the quark fields, Zi/2 = ("F^/f '^^'-'-^72-' lensing operator is 

2/' (22) = / dy^ [z2 ; y]8F'^'{y) [y; Z2] where, y^ =y^n^ +bj. For each Fock space state we have the phenomenological 
relation. 

We utilize this picture to calculate the gluonic pole matrix element using a soft approximation for the lensing 
function [20, 22] which we then convolute with the parameterizations of impact parameter GPD <§ [23, 24]. 
First we model the target remnant in terms of the sum of spectators the quark correlation function 
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2(27r)3(l-x)P+^ ^ 

which is expressed in terms of the matrix element W^'^ {P,k\a) = {P — k,a,5\ [°°n; 0]"^q'f (0) \P) where a and 8 
represent the helicity and color of the intermediate states (see Fig. 1). The FSIs- generated by the gauge link - are 
described by a non-perturbative amputated scattering amplitude {M)^g with j3, a (7, 5) color indices of incoming 
eikonal quark and out going spectator remnant. In momentum space W is given in terms of M 

AW«PfP,,, / ((P ,2^ [{F-4 + m,HP,SUM)fp{q,P~k) 

where i{n-{P — k — q) + iO) ^ ' represents the eikonal propagator and AW — W — W", where denotes the contribution 
without final-state interactions. Tracing Eq. 5 with 7+ and weighting, and integrating with respect to kj- yields the first 
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FIGURE 2. Left:The eikonal phase X^^{\zt\) vs. \zt\ with input from Dyson-Schwinger equations at scales Aqcd = 
OGeV, 0.2GeV, O.SGeV, 0.7GeV [20]. Right: The function C[f ] of Eq. (13) as a function of the eilconal amplitude |. We compare 
the numerical result computed by means of Eq. (10) up to the order « = 8 with the analytical result for the SU (2) color case. The 
numerical and analytical result agree up to ^ ~ 2. For SU (3), we compare the numerical results for the orders w = 7, 8. The results 
are reliable for j ~ 1.5. 



moment of the Sivers function [22] 



where the lensing function /' can be expressed in terms of the real and imaginary part of the scattering amphtude 
M [20]. 

We use functional methods to incorporate the color degrees of freedom for soft gauge boson coupling to highly 
energetic particles on the light-cone in the calculation of M, which is given by the expression 



aS , 



2{l~x)P+ -^^T^ 



(8) 



In Eq. (8) the — 1 dimensional integrals over the color parameters results from auxiliary fields a"{s) and u"{s) 
that were introduced in the functional formalism of Ref. [25] in order to decouple the gluon fields from the color 
matrices. The eikonal phase xdzrl) in Eq. (8) represents the amount of soft gluon exchanges that are summed up into 
an exponential form, and is given in terms of the gluon propagator 



Xi\zT\)=8 J daj d^n^n^^yiz + an-^n). (9) 

^ denotes the gluon propagator, and g the strong coupling. In this form the four-vector v is related to the complemen- 
tary light cone vector n, v = — ((1 — x)P^ / ms)n, with n-h=\ and = 0. We evaluate the color integrals by deriving 
a power series representation for the expression in brackets in Eq. (8), the color function f^^^ {%) 



{iXf 



f-p(x) - E E - E E(^"'-f""f"''""'-^"''""")a/3 ' (10) 

11=1 y"-> fl,=i fl„=i p„ 

where F„ represents the sum over all permutations of the set {!,...,«}. If we had a direct ladder where gluons 
were not allowed to cross we would have only factors {t"' ...t""t"" ...t"')ccp = C'p5ap with Cp = — l/2Nc, and 
we could work in an Abelian theory with an effective replacement a CfCCs for the fine-structure constant. Since 
we allow generalized ladders with crossed gluons we have to sum over all permutations in (10), and the sim- 
ple replacement is not possible. In a large A'^ expansion the crossed gluons diagrams would be suppressed such 
that the direct ladder represents the leading order in 1 /Nc- In an Abelian theory, the generating matrices t re- 
duce to identity and since we have n! permutations of the set {!,...,«}, we recover the well-known Abelian result. 
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FIGURE 3. Left: The lensing function /'(xJr) from Eq. (12) for U(l), SU(2) and SU(3) for x = 0.2 at a scale Kqcd = 0.2GeV. 
For comparison we also plot the perturbative result of Ref . [26] including the eikonalized antiquark spectator with an arbitrary value 
for the coupling, a — 0.3. Right: The first moment of the Sivers function versus x using various models for the GPD E. 



f^'-^Kx) = Lr=i ('Z)"/"' = e'^ - 1. For A^^^. = 2,t" = a"/! and we can calculate the integral analytically. We obtain, 
fftp^^ (x/4) = 5aj3 (cos;^/4 - ;t/4sin;^/4 — 1) + 18^1^ (2sin;^/4 + ;^/4cos;^/4). We also calculate numerically the 
lowest coefficients in the power series (10), and they agree with the coefficients in an expansion in X of the analytical 
result. This serves as a check of both numerical and analytical approaches. For A'^ — 3, due to difficulty of integrating 
over the Haar measure we use the power series (10) to obtain the approximative color function which is valid when 
a = ;)f/4 is small. 



- c^a — C(,a 



•)' = 5„p(cifl-C3fl^Vc5fl^-C7a^ + ...), (11) 



with the numerical values ci = 5.333, C2 = 6.222, C3 = 3.951, C4 = 1.934, cj = 0.680, ce = 0.198, cy = 0.047, 
eg ~ 0.00967. Transforming to coordinate space we can express the lensing function directly in terms of the real and 
imaginary part of the color function / which is itself a function of the eikonal phase % Eq. (9). This results in a lensing 
function of the form 

'X' 



^ ' ^' 2Nc \bT\ 4 



where we define the function. 
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(Ti-3 [/])'(! 
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(3[./])' (f ) m]) (f )] - ^Tr [(3[/]) (f ) imy (f 



(12) 



(13) 



and x' denotes the first derivative with respect to |z7-|, and (3[/])' and are the first derivatives of the real and 

imaginary parts of the color function /. 

In Fig. 3 the function C[j] is plotted versus j for various approximations. While the convergence of the power 
series seems to be better for SU(2) than in the SU{3) case where the numerical result calculated with eight coefficients 
agrees with the analytical result up to |- 2, we can trust the numerical result computed with eight coefficients up to 
f ~ 1.5for5t/(3). 

In order to numerically estimate the lensing function and in turn the Sivers function we utilize the infrared behavior 
of the gluon and the running coupling in the non-perturbative regime where we infer that the soft gluon transverse 
momentum defines the scale at which the coupling is evaluated. These two quantities have been extensively studied 
in the infrared limit in the Dyson-Schwinger framework [27] and in lattice QCD [28]. We use calculations of these 
quantities from Dyson-Schwinger equations [27] where both and f^^' are defined in the infrared limit (details can 
be found in a forthcoming publication). This determines the eikonal phase and thus the lensing functions (12) for a 
t/(l), SU{2) andSU{3) color function. We plot the results in Fig. 3 for a color function for f/(l), SU{2), SU{3). While 
we observe that all lensing functions are attractive and fall off at large transverse distances, they are very different in 
size at small distances. 

Using the eikonal model for the lensing function together with different models for the GPD E we present the 
first moment of the Sivers function from Eq. (4). In Fig. (3) we compare the diquark model result [22], to two 



phenomenological models 1 24, 23] . We find that the lensing function grows with Nc which in turn predicts the growth 
of the Sivers function. This is consistent with the large Nc QCD scahng behavior of Sivers functions [29] and further 
we uncover the deviations for finite Nc. Further, the sign and size of the Sivers function is consistent with present 
extractions from data [30]. Within this framework of FSIs we have performed quantitative analysis of approximate 
relations between T-odd TMDs and GPD which goes beyond the discussion of overall signs. 
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